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Generalization of the DCX Reaction Amplitudes 
for Transitions to an Arbitrary Final State 


Uogólnienie amplitud reakcji DCX dla przejść do dowolnego stanu końcowego 


1. INTRODUCTION 


In a series of papers [1-4] (hereafter referred to as KF) it was shown 
that the double charge exchange reaction (DCX) with pions can be descri- 
bed within the proton-neutron quasiparticle random phase approximation 
(QRPA). Both transitions, to the ground state and to the double isobaric 
analogue state of the final (A, Z + 2) nucleus were studied and the expres- 
sions for the total amplitudes in these cases were found. The question arises 
as to whether such a result can be reached if we take into account the trans- 
ition to any arbitrary excited state of the final nucleus. We shall investigate 
this problem in this paper. 


2. KF MODEL REVISITED 


In order to generalize the KF formulae, it is instructive to review briefly 
main goals of the previous results. We have considered the total transition 
amplitude F(k, k’) for the DCX reaction 


at + X(A,Z) — x7 + Y*(A,Z +2) (1) 


defined in such a way, that the differential cross section is equal to 


do 


AS = ąz Füe kf. (2) 
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Here k and k’ are momenta of in- and outcoming pions, respectively. We 
applied to reaction (1) the sequential model in which an incident pion r+ 
changes in a first step into a pion 7° by changing a neutron into proton and 
producing the intermediate nucleus (A, Z+1) from the target nucleus (A, Z). 
During this process, the intermediate nucleus is excited to all possible 
states. In our model, these excited states are described in the proton- 
particle neutron-hole QRPA formalism. We did not involve the effects of 
a nucleus penetration by the neutral pion explicitly because they are higher 
order terms in the linear response function. Next, a pion x~ is produced by 
changing another neutron into a second proton and leaves the final nucleus. 

The total amplitude F(k, k’) is a sum over all excited states of the 
intermediate nucleus. We restrict our considerations to nucleon excitations 
only in the low pion energy domain. In principle there are no fundamental 
difficulties to include delta-isobar excitations, but because of the large 
delta-nucleon mass difference such terms are of no importance for pion 
kinetic energies lower than approximately 150 MeV. So, the full scattering 
amplitude is 


F(k,k') zi > [Fos(k, k’)] , (3) 
{aJ} 


where we sum over all intermediate states (a, J} as stated above. To de- 
scribe the amplitude (3) in a more detailed way we introduced two pheno- 
menological transition Hamiltonians reflected two different mechanisms for 
the DCX reaction. One of these Hamiltonians is generated from the non- 
relativistic reduction of the well-known pseudoscalar coupling Lagrangian. 
In this way one can obtain so-called p-wave effective pion-nucleus Hamilto- 
nian (for details see [2]) 


hp(q) = -v2 i JY e , (4) 


Mr pn JM 


Another mechanism is connected with an exchange of a composite boson 
between nucleons. The mrN interaction playing a crucial role in such a 
process gives us the phenomenological s-wave rN effective Hamiltonian i in 
the Koltun-Reitan’s form [2] 


A 
hula) = -4r <V 2w, jp Rum, (5) 
c pn 
In the above equations the QRPA transition operator a is defined by the 


formula 
RIM = wv„Ct(pnJM) + tptnC (pnJ M)+ 
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+ UptnD(pnJM)— vjv,D'(pnJM) . (6) 


Here w, = (q? + m?)! is the incident or outgoing pion energy, j= 27 F 1 
and u, v's are the occupation amplitudes appearing in the transformation 
from the particle operators to the quasiparticle operators. The pair-creation 
(annihilation) operators C'(C) and additional Dt(D) operators are defined 
in the usual way (see e.g. [2, 5]). The function FM (q) can be written as 
follows 


FIM(q)= XO (-1)*77%*(jpmpjn — m4|JM)x 


Mp,Mn 


x || Przygoda aee] (7) 


The proton (neutron) space-spin wave functions pp(n)(x,€) are calculated 
with some single particle potential (e.g. the harmonic oscillator or Woods- 
-Saxon potential). After some algebra one can fina a more handy expression 


FIM (q) = Varv6Yjy(Qa)Gyn(a) . (8) 


where Yyy is the spherical harmonic depending on the solid angle 2, and 
G7, (a) is the nuclear form-factor, for which the final expression is different 
for different choices of the radial dependence of the nucleon wave functions. 
More details on the explicit form of Gz, can be found in [2]. 

Using perturbation theory we could express the transition amplitude 
Fj(k, k') as follows: 


x” (k'|f|mJ M)(mJ M|T|i, r+ (k)) 


Fj(k,k') = |) if P" GE (9) 


where T is the transition operator consisting of both s- and p- parts of the 
pion-nucleus Hamiltonian (eqs. (4) and (5)). |i,r*(k)) is the ground state 
of the target (A, Z) nucleus with the incoming pion t+ of momentum k. 
( E; + wk) is the initial energy of the nucleus-pion system and E;,'s are the 
QRPA excitation energies of the intermediate nucleus. In eq. (9) we sum 
over all intermediate excited states. In point of fact, the question of this 
summation is more complicated because of two different mathematically but 
equivalent physically sets of the intermediate states calculated relative to the 
target nucleus (A, Z) and the final nucleus (A, Z + 2), respectively. We do 
not intend to discuss this problem further in the paper. 

In the approach developed earlier, the state |f,r”(k')) considered was 
either the ground state of the final nucleus or the double isobaric analogue 
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state in the same nucleus. In principle there are no fundamental difficulties 
to treat the final state as an arbitrary excited state of the (A, Z+2) nucleus. 
A need for such consistent formulae is urgent all the more, that the DCX data 
exist for the reaction 5$Fe(rt,m”)5*Ni in a full extension. Cross sections 
and angular distributions for the double isobaric analogue transition, non- 
analogue ground state transition and for some transitions to the low excited 
states (07 ,OŻ ,27 ) were measured recently [6,7]. 


3. GENERALIZED DCX TRANSITION FORMULAE 


The state |f) in the expression for the total amplitude (9) can be 
specified as the ground state (GS) of the final nucleus or one of the excited 
states among which the double isobaric analogue state (DIAS) is mostly 
discussed in literature. In KF we derived formulae for the DCX transition 
to the ground state [2, 4] and to the DIAS [1]. Below we shall show how 
to generalize such formulae in the case when the final nucleus is excited 
during process (1) to the one of its excited states. We are working within 
QRPA formalism to describe both sets of states: the excited states of the 
intermediate nucleus and the excited states of the final nucleus. Thus, 
we have to define two types of the angular-momentum-coupled phonon 
operators from the usual RPA ansatz [8]: 

(19) for the excited states of the intermediate nucleus 


UWM = p= [XGnyC'(pnJM) = Y n)sĆ(pnJM)] : (10) 
pn 


(29) for the excited states of the final nucleus 


Qu! = Y, Eos App I M) - Kępy m AGE I M)] + 
pp’ 


+E [s Bin’ IM) = Vian gaeB(nn'TM)| „| (11) 
where pair creation (annihilation) operators for the proton-neutron system 
C'(C), proton-proton system At(A) and neutron-neutron system Bt(B) 
are defined in a usual way. X's and Y's are forward- and backwardgoing 
QRPA amplitudes of the mth J state of the odd-odd intermediate nucleus 
calculated relative to the final nucleus ground state with proton-particle 
neutron-hole excitations (the charge-changing mode). Similarly, ¥ (ear) and 
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Vac’) with (aa') standing for proton-proton or neutron-neutron indices are 
the QRPA amplitudes of the excited states in the final even-even nucleus 
calculated within the charge-non-changing mode. According to eq. (9) we 
need to calculate two matrix elements 


(mJM|T i, GS; x* (k)) = (i, G5|QTyT li, GS; x* (k)) (12) 
and 
(fur IM; x7 (kITImJM) = (,GS; 7 (K))O s TQ GS) , (13) 


where JM is the angular momentum and its third component of the 
intermediate state and JM is the angular momentum and its third 
component for the excited state of the final nucleus. In real applications 
of the formalism we use two different phonon operators of type (10), 
because in the charge-changing mode one can calculate the excited states 
of the intermediate nucleus in two mathematically different but physically 
equivalent ways. From a formal point of view the only difference is that we 
should distinguish the RPA amplitudes X and Y's by adding a bar in the 
case we calculated them starting from the ground state of the final nucleus 
or leaving a bar in the case we calculated the amplitudes relative to the 
ground state of the initial nucleus. In this point a more detailed discussion 
of such difference is not needed and does not influence on generality of final 
results. 

To find formulae for the matrix elements (12) and (13) we should 
remember that the QRPA phonon annihilation operators Q zy and Q^, 
give zero operating on the RPA vacua. We also use the quasiboson 
approximation [5] and finally we are taking the BCS solutions as the RPA 
vacua. With the above approximations one finds after little algebra: 

(19) for the s-wave transition operator h, 


(RPA| [R99 Ct(p'n'JM)] |RPA) = tpunbpyrównójoómo >» — (14) 
(RPA| [Rpa, C(p/n'JM)] IRPA) = (-1) * * upon bzpónnójoómo „ (15) 
(RPA| [[4(/7Mm),ff,.] , Ct (pn M)] IRPA) = 

= uyrusn( 1r In 1 ją (Gypróprpónn" Sin jp t 


+(-1) dypnOppóznója;» )BJTÓMM > (16) 
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and 


(RPA| |[At(pp' 7 M), RA,.] ,C(piJ M)] [RPA) = 
= vpnvan(—1)# rH jn (bpp Spry bin Sing + 


HL) dygróypównnó;1;,) 275MM » (17) 


Similar equations with the neutron-particle neutron-hole excitations can be 
written immediately by exchanging proton indices on neutron indices. 
(2°) for the p-wave QRPA transition operator RZ! i 


(RPA| Rzut" , C'(pnJ M)] [RPA) = 


= (71) *M opnttgn bpp San "MMe 3 (18) 


(RPA| [774 ,C(pnJ M)| (RPA) = 


= (- 1)” ar Up Vn" Sop" San Oar" Ons age ; ( 19) 


(RPA| ([A(pp'T.M), RM") ,C' pas M)] (RPA) = 


= yay (- Lyr JPG M Im i Tw 5 [p 
| P 


* * J , 
"e dt. Je L b. aluit 
(7-1) y iy p'p""pp^R 


: (20) 


and 
(RPA| ||At(pp' JM), R^] ,C(pad M)| [RPA) = 


= vpron ( -t Fr gues [s ja T ü barb opns T+ 
Pp 


ar (-1)7 je p re ) d L (21) 
In expressions (14) — (20) the symbol [. , ] means the commutator and 
daja! = ÓnanaÓj ju Ślal„ Above we only quoted the terms with non-zero 
contributions within the quasiboson approximation. 

With an aid of formulae (14) - (21) one finds the needed matrix element 
(13) containing two parts for h, and h, , respectively: 
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(feJ M;z”(k')|h,|mJM) = 4x A v2 0X6] 36M MX 


- * sa m a m iT . 
x y: (-1yrt^ [Er m cm Vipp) Y pns Yva) bp 4 pat 
PSP! pn 


+ (-1)’ (35, AA = Pips fga 3% 9") öpp] 9 (22) 


UJ M; T (k')|h,|mJ M) = iL Vas Tix 


zen JII MI" M" My cay Yj |f Jo ah 
PSP! pad" M" Jp 


" -= > FTR = > Las j ] J 
X6ygGj; (k^) (As X nap E Vopr) s¥ mys Bpin ) + m pA d ) X 
J" po Pm = = V * -— 

x (71) 5567, (/) (A, nirin — psi urs m)| . (23) 
In all above expressions the quantities with bar correspond to the states 
generated from the daughter (A, Z + 2) nucleus. Both formulae (21) and 
(22) are quoted for the proton-particle proton-hole part of the total matrix 
elements only. The discarded neutron-neütron part of eqs. (22) and (23) 
may be easily constructed by exchanging p, p' indices upon n,n’ indices. 

Analogous equations can be written for the states generated relative to 
the initial (target) (A, Z) nucleus using the RPA a; plitudes without bar. 
Thus, the matrix element (12) can be written as foliows [2]: 


(mJ M|hili, G.S. r*(k)) = v2 (a>) TRU 


= 
and 
(mJ”"M|h,|i,g.s., x *(k)) = i V12 L Var (-1)7 Yžu()x 
x |2/G>(k) (X7jtpta i Xm) (25) 
.pn 
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Combining eqs. (22), (23), (24) and (25) with the expression for the total 
amplitude (9) one can find the final expressions for the s- and p- wave parts 
of the DCX amplitude in a case of the transition to an arbitrary excited 
state |J M) of the final nucleus: 


Y (mJ|m'J) r 
J 
m,m’ Ei + wę — Satin 


4 > 42 
Fi(k,k') = —ÓJ7 (405) WĘWĘ: 


x (ve E > (Eor X Gems etn - i a (sn) pln 
PhupSp' 


+(-1)7 (Xr X (nya fign ^ ‘op. (yyy da) Gppózya | x 
x [V2 830925 Son (Xen rupon — ¥(™%) tpun )| ) , (26) 
pn 


F5(k,k’) = — (2) Ka a (mlm) x 


E),+E7, 
Mz m,m’ E; + ty — mt 


x 5) II-T HM Ym) Yd] 
JH M xk 


Vi2 Y^ (-iyrt?x 


PSp’ pn 
dtd - pay bo. eating 
i ii I" jp } Sys Gon (K) (Zip X conyatipiin — Vege s¥ guys pn) + 


Tay 
{HN DCC 


A (Fors Xos ewa Yis | á 


x [vi E Cim (k) (XGpuyattp Pe — rw)! ZU 
pn 


In eqs. (26) and (27) we involved a summation over both sets of the 
intermediate states |mJ M) and |m'J M) generated independently from the 
parent (A, Z) nucleus and from the daughter (A, Z 4-2) nucleus, respectively. 
They fulfil the orthogonality relation only approximately, so we are forced 
to introduce overlaps (mJ M|m'J M) into the final expressions for the 
amplitude (9). In addition, the procedure introduces an uncertainty in the 
energy denominator and we decided to adopt an average value of the QRPA 
energies EJ and EJ, instead of one of them. In general, this uncertainty is 
small however, compared to the mean value of the denominator. 
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4. SUMMARY AND CONCLUSIONS 


In this paper we have generalized the formulae for the DCX transitions 
to any excited state of the final nucleus. As in the previous papers we used 
the proton-neutron, proton-proton and neutron-neutron QRPA to describe 
consistently both excitation modes (charge no — and — changing ones) 
appearing in the reaction model. The expressions obtained are in agreement 
with the previously derived result for the final double isobaric analogue [1] 
state if one assumes the angular momentum of the excited state in eqs. (24) 
and (25) equals only to zero. 

The generalized result will be applied to the iron and nickel isotopes in 
which complete data for three types of transitions exist [6]. Because different 
possible routes influence on the DCX cross section [2,9,10] depending on 
the type of final state (GS, DIAS or another excited), one can expect to 
have in such calculations a very promising way to study various short-range 
correlation effects between nucleons in the nucleus. Such a project is being 
developed now and results will be published elsewhere. 

This work was supported by the grant from the Polish Ministry of 
National Education. 2 
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STRESZCZENIE 


W pracy uogólniono wyrażenia na przekrój czynny w reakcjach podwójnej wymiany 
ładunku z pionami. Rozważono przypadek, w którym jądro atomowe zostało wzbudzone 
w trakcie procesu do któregoś z jego stanów energetycznych. Przytoczono formuły na 
różniczkowe przekroje czynne jako funkcje kąta rozpraszania 0 dla rozpraszania w fali si 
w fali p. 


